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ABSTRACT 









A sea floor benchmark experiment was conducted in an 
area about 16 nmi west of Pt. Lobos, California (36°30'N x 
122°17' W) during 18-22 May 1985. 

Two baseline-crossing methods were used to determine the 
relative positions of acoustic bottom transponders . The 
method of least-squares adjustment was used to analyze the 
data. Relative position determination of the transponder 
array is discussed and recommendations are made for futher 
improvement. The advantage of these methods is their 
simplicity. Their disadvantage is the relatively large 
amount of ship time they require to achieve acceptable accu- 
rac i es . 

Transponder arrays such as the one deployed can be used 
for solving many types of problems in sea floor engineering, 
which will be of increasing importance in the future. 
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I. I MTBODOCT ION 



A. BACKGfiOOND 

In today's world, the ocean scientist and engineer need 
increased accuracy for marine positioning to accomplish 
their work- The establishment of a system for the precise 
determination of positions on the deep seafloor is one of 
the most fundamental challenges of working in the open sea. 

A seafloor benchmark positioning system is proposed 
which can aetermine accurately the location of- objects on 
the sea floor (Saiena, 1974). Such precisely located benca- 
marks can later be used in turn to delineate offshore prop- 
erty lines and national boundaries- They are also useful in 
solving problems in seafloor engineering, in plate tectonic 
studies and in connection with bore-hole reentry associated 
with off-shore oil recovery. 

The objective of this study is to demonstrate the feasi- 
bility of establishing such seafloor benchmarks. 

B. DESCEIPTIOM 

Early experiments associated witn seafloor acoustic 
transponder arrays were carried out by Hart (1967), Haehnle 
(19o7), Fubara and Mourad (1972), and McKeown (1975), who 
used a baseline crossing method to solve for the distances 
between transponders. "Baseline" in this thesis means the 
line between the projections of two bottom-mounted acoustic 
transponders onto the sea surface- 

Ihe baseline accuracies obtained by Hart (1967), acKeown 
(1975), and Fubara and Mourad (1972) were ±15.7 m for base- 
lines of 5509 m, ±3. 1 to ±4.2 m over baselines of b373 to 
7219 m and ±15.5 m over baselines of 9364 m, respectively. 



The first two of the above authors did not use a least 
squares method for data analysis- 

There are three parts to the determination of the abso- 
lute position of an ocean-bottom transponder array using 
satellites- They are (Figure 1-1): 

1- To determine the geodetic position of the ship's 
receiving antenna from the Global Positioning System 
(GPS) ; 

2- To convert the geodetic position of the antenna to 
the ship- mounted acoustic transducer; and 

3- To determine the geodetic position of the acoustic 
transponder array on the sea floor using a snip- 
mounted acoustic transducer- 

This thesis is concerned with the relative position of a 
bottom- mounted transponder array with regard to the ship- 

C. DESCBIPTION OF THE SYSTEM 

Acoustic transponders manufactured by Oceano Instruments 
were used for this project. They were deployed by the R/V 
Acania. Figure 1.2 shows the configuration of the Oceano 
system used- ' Their principal features are: 

1- AM-121 : (Acoustic Module) 

The hydrophone was fixed amidships to the 
port side at a depth 2. 95 m- 

- Transmission frequencies : 8 to 16 kHz- 

- Reception frequencies : 8 to 16 kHz. 

2. HT-121 : (Recoverable Transponder) 

- Frequencies : 8 to 16 kHz- 

- Delay time : 15.0t0.1 ms- 

TT-201 : (The Accoustic Telecommand Module) 

has the following functions: 

- Code acoustic signal transmitted by either the 
acoustic module or UQC transducer. 

- Monitors range to desired transponder. 

1 1 
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Figure 1- 1 



The Configuration of Bottom Transponder Array 
Shown with Ship and GPS Satellite. 
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- Monitors reception of code by transponder. 

IM-100 : (Data Interface Module) 

- Accepts up to 16 data inputs. 

- Max. converted data format : 7 digits without 
sign or 6 digits plus sign. 

- Computer interrupt request possible through 
any input. 

- Each channel may be read separately. 

5. RM-201 ; (Eangemeter) 

- Transponder turnaround delay offset : 0 to 99.99 ms. 

- Mean velocity of sound input : 0 to 9999.9 units per 
second. 

- Standard reception filters : 8 to 16 kHz {0.5 kHz step) 

- Bandwidth : 200 Hz. 

Critical to the experiment are the following: 

1. The transducer (a part of the AM-121) was 2.95 m 
below the surface, and the transmission and reception 
frequencies were 8 to 16 JcHz. 

2. The system delay on the EM-201 was set at 15 millise- 
conds, i.e. , the time lapse between reception and 
transmission of a siqnal. 

3. Velocity on the RM-201 was set to 1480 m/s. 

4. Interrogation period was set at 10 s. 

5. Interrogation frequency (kHz) was set 15 kHz. 

The numerical constants of (3) and (4) were used for the 
survey. Details of the transponders are given oy Oceano 
Instruments (1984). 
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Figure 1.2 



Configuration of The Oceano Systems Used 



II. THE OPERAIIQBAX PROC EDURE FOR I MPL ANTAT ION OF A 

TRANS POND ER ARRAY 



A. SITE SELECTION 

The main criterion for the site of the experiment was 
that it be relatively flat over an area approximately 2 km x 
2 km. Since there were no existing charts detailed enough, 
a reconnaissance survey was perfomed by the R/V ACANIA on 10 
and 11 April 1985. Some 10 north-south sounding lines and 
1-cross line were run with 1-nmi spacings and 5-min position 
fixes. 

The area chosen for the experiment (Figure 2.1) is 
approximately 24 nmi west of Pt. Lobos, California, in water 
approximately 1500 to 2000 m deep. It is the area from 
36°29.5' to 36<>31'N and 122016 * to 1220 17.8'H. Mini -Ranger 
stations (Figure 2.2) were set up at Point Sur 
(36018* 20.279“ N X 121o 53* 56. 179"W) and Carmel 
(36033* 49. 176'* N X 1 2 1® 53 * 4 8. 3 58**H) for additional ship 
position control. 

B. IMPLANTATION PROCEDORE 

The transponders were deployed at pre-selected sites 
using a buoy-first, anchor- last, technique. The anchors 
used were railroad wheels which weigh about 290 kg in water. 
The mooring cables between the anchors and the transponders 
and between the transponders and the buoys are each 8-m 
long. The transponders weigh about 46 lbs in water, and 
each buoy has 55 lbs positive floatation. Three 17-in 
Benthos glass sphere floatation buoys were attached to each 
transponder by means of 5/32-in stainless cables 8-m long. 
Altnough 2 spheres would have provided sufficent buoyancy 
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Figure 2- 1 The Project Area. 




Figure 2.2 The Project Area with Shoreline. 
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for recovery, a third was used to provide additional tension 
(about 50 kg total) to reduce the motion of the transpon- 
ders. The mooring system is shown in Figure 2.3. 




Schematic of Bottom Transponder Mooring. 
1 9 



Figure 2.3 



III. DITEBMIMTI£H of ha h mon ic mean soond velocity 



A. INIBODUCTIOH 

Historically, echo soundings have been made by assuming 
a constant, approximately average, value for the velocity of 
sound throughout the water column, usually 
4800 ft/s (=1463 m/s) or 1500 m/s, and then correcting for 
variations of this assumed value for the actual water 
column. Echo sounders are time-measuring devices. The 
sounder's acoustic transducer emits a sound puise which is 
reflected upwards by the bo ttom- mounted acoustic transponder 
and received back at the surface transducer. The time of 
travel of the sound pulse is divided by two, and this value 
is multiplied by the assumed harmonic mean speed of sound in 
sea water, thus, giving the distance according to the 
expression D = V t . The echo sounder makes this transfor- 
mation electronically or mechanically within the device 
itself and displays the aistance in the water between trans- 
ducer and transponders. This distance is not equal to the 
true distance, since the assumed harmonic mean sound speed 
generally does not equal the true mean speed for a partic- 
urar depth. To determine the true harmonic mean sound speed 
one must know the sound speed throughout the water column 
and apply it to solve for the harmonic mean. For a distance 
measurement to be accurate, precise measurements of time and 
the harmonic mean speed of sound in the water column are 
mandatory. 

B. THE VELOCITY OF SOOND IN THE OCEAN 

The velocity of sound in sea water is a function of 
temperature, pressure, and salinity: V = V(T,S,P). Because 
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tne sound velocity is not constant with position, sound rays 
are refracted according to Snell’s law. 

We have followed the standard practice of the U- S. 
National Ocean Service (NOS) to use Wilson’s (I960) equation 
for the speed of sound as a function of salinity, tempera- 
ture and pressure (Umbach, 1976) i 



V=1449 . 14+V^+Vp+V2+Vg,pp 

V^=4 . 5721T-4 . 4532X10“^T^-2 . 6045X10"^T^ 

+7.9851X10"®T^ 

Vp=l . 60272X10“^P+1 . 0268X10“^P^ 

-(-3. 5216X10"^P^-3 . 3603X10"^“P^ 

1 

Vg=1.39799CS-35)+1.69202X10”^(S-35)^ 

VsTp=( S-35 ) ( -1 . 1244X10'^T+7 . 7711X10" V 

+ 7. 7016X10"^P-1 . 2943Xlo’^P^+3 . 1580X10 
+1. 5790X10"^PT^)+P(-1 . 8607X10"^T 
+7 . 4812X10"^T^+4 . 5283X10"®T^ ) 

+P^C-2 . 5294X10"’^T+1 . 8563X10"®T^ ) 

+P^C“1 . 9646X10"^°T) . 



(3.4) 

-8pT 

(3.5) 



In these equations the absolute pressure, P, is expressed in 
kg/cm2, temperature, T, is in °C, salinity, S, is in g/kg, 
and sound velocity, V, is in m/s. For these relations for V 
the standard deviation from the mean is 0-30 m/s for all 
data obtained in the ranges -4 < T < 30 °C, 1 < P < 1000 
Jcg/cm2, and 0 < S < 37 g/kg (Wilson, 1960, p. 1357). 

C. HABMONIC MEAN SOOHD VELOCITY 

The slant range from a ship to a transponder is deter- 
mined by measuring the transit time, t, of an acoustic pulse 
and converting it to distance using an appropriate value for 
sound speed. If the mean sound velocity, V, in the water 
column is known, the distance, L, from transducer to tran- 
sponder can be computed by: 

L = ( t/2 ) V (3.6) 

For our small project area the distances between transpon- 
ders are less tnan 2 nmi, and it is assumed that there is no 
horizontal sound speed variation, and only corrections for 
vertical sound speed variation are considered. For a 
one-way travel time , T, through a water column of depth, Z, 
the mean sound velocity is: 

V = Z / T (3.7) 

The meaning of '•mean sound velocity" is neither "mean 
velocity from the surface to the stated depth'* nor the "mean 
value for the velocity of sound through the vertical water 
column". According to Equation (3.7) the mean velocity is 
obtained from toe true depth divided by the one way travel 
time- 
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Figure 3-1 A Layered Ocean Bhere ?;= (z) Is the Sound 

Speed in the xth Layer and z Is Depth- 

Consider nov a layered ocean having a different sound 
speed in each layer. Figure 3,1 shows a series of finite 
layers of thicKness az. , each with an associated sound 
speed, . The time interval, AT., for the sound wave to 
pass vertically through the i layer is : 

AT- = Az^ / 7^ (3-8) 

Summing tne time intervals for all n layers in a vertical 
water column from surface transducer to bottom transponder , 

T = r at. = £ (AZ. /V^ ) (3.9) 

= f Uj 

For a continuous function, V = V(z), this becomes 
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T 




dz/V(z ) 



(3. 10) 



Equation (3.7) may now be written : 



V 



Z 



(3. 11) 



i. 



Uz 



V(z) 



This is the integral form of the harmonic mean- Using 
Simpson’s Rule, Equation (3.10) can be rewritten as : 



where d = D/n (magnitude of depth increments) 

D = bottom depth at transponder 
n = even number of depth increments between 
transducer and transponder 
V; = measured sound speed at depth j. 

Since the harmooic mean sound velocity is obtained by 
Equation (3.11), the slant range from transducer to acoustic 
transponder can be obtained using Equations (3. 10) and 
(3. 11): 

L = T V 

CTD casts to 140 0 m were made on 18 and 20 May 1 985 in 
toe project area just prior to and during our main experi- 
ment. Sound velocities and harmonic means based on the 
sound velocity profiles derived from the CTD casts are 
plotted in Figures 3.2 and 3- 3 . Altnough the profiles 
shown in these figures are for different days and different 
places within the project area, they show very little varia- 
tion in harmonic mean sound velocity, and thus support our 
previous assumption that the sea water is horizontally homo- 
geneous to an acceptable degree within the study area- 



T = ( 1/V„ +4/V, + 2 /V 2 +4/V^ + ••• +I/V 2 ) d/3 



(3. 12) 
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Figure 3.2 Harmonic Mean and Sound Velocity Profile 1. 
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Harmonic Mean and Sound Velocity Profile 2. 
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Figure 3- 4 Hacnonic Mean Cooparison. 
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Actually, the difference of these two harmonic mean veloci- 
ties is approximately 0-07 m/s near the bottom as is shown 
in Figure 3-4. The mean of these two harmonic mean sound 
velocity profiles is used in the range calculations of the 
following chapter, 

D, ACOOSTIC WAVE EEFBACTION 

To simplify our acoustic analysis it is desirable to 
assume that in the study area sound follows straight paths 
from the surface to the ocean bottom- To test this approxi- 
mation the differences in path lengths from the surface to 
the bottom were calculated for varies depressi-on angles at 
the surface- 

profiles such as those of Figures 3.2 and 3-3 are 
usually simplified for analysis by separation into an appro- 
priate number of segments each having an approximately 
constant gradient (Figure 3.5). If the velocity of sound 
changes linearly with depth, sound rays can be shown to have 
a constant radius of curvature (Uricic, 1983, p. 124). 
Horizontal distance, d, depth, z, radius of curvature, R, 
segments of arc, s, and depression angle. A, are shown sche- 
matically in Figure 3.6. The path of a ray when the speed 
of sound varies witn depth can be calculated by application 
of Snell's law: 

cosAj / C| = cosA^ ^ ~ *** ~ cosA^/ (3.13) 

wnere C*= sound speed 

Kl= depression angle 

For an initial depression angle of 30^^ the profile of 
Figure 3.5 gives the results shown in Table I. 

The raaius of curvature can be obtained in the following 
manner (Kinsler, 1982, pp399-402). Let 
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Figure 3-5 Simplified Sound ?elocity Profile- 
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Figure 3.6 Schematic Diagram of Sound Ray Trace. 
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TABLE 


I 








Depression Angle at 


Different 


Layers 




1 


Sound 1 


Depth 


1 


D epression 




Points 1 


Velocity | 


(in) 


1 


angle (deg.) 




1. 1 


1493.0 1 


0- 


1 


30. 0 ‘ 




2. 1 


1483.5 1 


176.4 


1 


30.63 




3. 1 


1484.04 1 


266. 6 


1 


30.59 




4. 1 


1481.2 1 


526.4 


1 


30. 78 




3. 1 


1482.46 I 


1097.2 


1 


30-69 




6. 1 


1484.34 1 


1400. 0 


1 


30.57 


5 










• 




= < - 


- 






(3. 14) 


and 












B- 

L 


= - (C. / 


g. ) cosA- 

L ^ 






(3. 15) 


where 














= Radius 


of curvature in 


the i ^ 


depth increment 



C; = Sound speed in the depth region 

= Depression angle 
g. = Speed gradient 

The change in range Ad- and depth az£ are (Kinsler, 1982) : 
Ad|^ = -R ^ ( sinA^ - sinA^^^j ) (3.16) 



= -E^ ( cosA^^j - cosA^ ) 



A z- 

L 



(3. 17) 



so the chord can be obtained by ttfo components of change in 
range and depth; that is, 

D. = ( ^d,2 + Az.2 )i/2 (3.18) 

t c 

The angle, 0- , is 

9; = 2sin-i ( D; /2R; ) (3.19) 

I- L ^ 

so the arc, s^, is: 

s • = R- 9* (3.29) 

l- L c - 

where 9^ is the angle in radians. 

The range and depth increments give the straight line 

L = [ (lAdi ) 2+ Uaz- ) 2]1/2 (3.21) 

The ray paths can be obtained by summing arcs: 

S = £si (3.22) 

The difference due to refraction between a sound ray's 
actual path and a straight line is thus: 

E = S - L (3. 23) 

where E = the error between two lines, 

S = the distance along the ray path, and 
L = the straight-line distance between the 

ship-mounted transducer and bottom transponder. 
Ray-trace computations were made for various initial 
depression angles (at the surface) down to a depth of 1400 m 
(Taoles II and III)^Since E<1 cm for a 30® surface 
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TABLE II 

The Result of Ray Trace Computation 





For Initial 


Depression Angle 


= 30 Degrees 


Depth Gradient 


Radius 


Horiz. 


Depth 


Arc 


Increment g 


1 

1 

1 

1 

1 

1 


Dist. (d 

.JjnJ 


) z 

Lnil... 


1 

{ 

1 

1 

1 

1 

1 


1. 


0. 0538 


32011. 36 


301. 72 


176.4 


349. 5023 


2. 


0. 006 0 


287966, 49 


152.47 


90,2 


177. 1543 


3. 


0. 0109 


157706..64 


437,85 


259.8 


509. 1278 


4, 


0,0022 


780984- 82 


960.03 


570, 8 


1116, 8993 


5. 


0. 0062 


277668. 87 


511-35 


302.8 


594. 2820 




Total : 2363.42 


1400.0 


2746. 9658 




The Distance 


’of Straight Line 


• 

• 


2746. 9562 




The Error 


• 


0. 


0096 mete 


rs 



depression angle and a 1400-m water depth, corresponding to 
a horizontal distance of 2363 m, which covers the project 
area, and since the depression angles encountered during the 
project were 30° or greater, 1 cm represents an upper limit 
to the error resulting from our assumption of linear ray 
paths from the surface to the bottom, and a correction for 
refraction is unnecessary. 
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TiBLE III 

The Error Coaparation at Different Angles 



1 



Depth : 1400 meters 



De pression 
Angle 
( degree) 


Horizantal 

Distance 

(in) 


Ray 

Trace 

(m) 


Straight 

Line 

(m) 


Error 

(cm) 


0 


13429-23 


13508.3850 


13502. 0064 


637.86 


10 


b637 .03 


6783.2984 


6783.0795 


21. 89 


20 


3647. 10 


3906.6140 


3906.5806 


3. 34 


30 


2363.42 


2 74 6- 96 58 


2746.9562 


0. 95 


60 


801 . 28 


1613,0859 


16 13.0853 


0. 06 


90 


0.0 


1400.0 


1400.0 


0. 0 
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IV. R ELATIVE POSITION PETEK MIN ATION OF THE TRANSPONDER ARRAY 

A. INTRODUCTION 

Two methods to determine relative positions are 
discussed in this chapter. They are basic baseline crossing 
and the modified baseline crossing method; the former is 
described by Hart (1967) and Haehnle (1967), the latter is 
described by Fubara and Mourad (1972). The main difference 
between the methods is that Hart and Haehnle reguire accu- 
rate transponder depths, but Fubara and Mourad- do not. 

The cloverleaf method used to determine transponder 
depths IS described by Haehnle (1967) is introduced in tne 
next section. 

B. CLOVERLEAF METHOD FOB TRANSPONDER DEPTH DETERMINATION 

An initial approximation to the depths of bottom-mounted 
acoustic transponders can be obtained from a knowledge of 
tne water depth at the point where the mooring is released 
from the ship. However, the unit may move laterally while 
descending, or bottom topography may have considerable 
relief not evident on the echo sounder record- Hence, a 
more accurate method of depth determination is required. 

Here the depth is defined as the distance between the hull- 
mounted transducer and the bottom-anchored transponder when 
the vessel is positioned exactly over the transponder- 

The ship need not be directly above the transponder but 
should be in close proximity to it. Although it is very 
difficult to cross over the transponder exactly with tne 
ship, the depth may be approximated as the minimum slant 
range when the transponder is crossed over many times using 
a '• c 1 overleaf *' method (Figure 4,1). Thus, if a transponder 
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a. The firs-t pass through the transponder location. 

I 

I 




b. Maneuver for obtaining four CPA's. 

i 

I 



Figure 4«1 Deteraination of Transponder Depth- 
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is at a depth of 1530 m aad the ship is dxspiaced 50 meters 
horizontally from it, tne depth determined by the cloverlea 
method would be in error by only 0-33 meters (Figure 4.2). 

The procedure used for running a cloverleaf maneuver is 
as follows (Haehnle, 1967) : 

1. Each time a transponder is dropped from the ship the 
ship's position is noted. 

The ship proceeds toward the transponder drop point 
and passes near it; A closest point of approach 
(CPA) is then determined and recorded- 
The ship then turns, crossing the first track perpen 
dicularly and passing near the first CPA. Once 
again, the CPA of this run is recorded- 
The above maneuver is repeated until at least four 
such CPA's have been obtained, each time attempting 
to pass through the previous CPA perpendicular to th 
track. 



2 . 



3. 
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Figure 4.3 Saiple of Depth Determination- 

5. Ihe depth of the transponder is deterained by plot- 
ting the slant ranges versus time for each of the 
four CPA runs. The minimum point of this curve is 
considered as the minimum slant range for that 
particular transponder (Figure 4.3). ' 

The R/V ACAilIA ran similar cloverleaf patterns on 21 May 
1985. She crossed near each transponder many times as 
closely as posstnle. Ihe minimum distances and corrected 
depths are listed in Table IV . 
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TABLE I? 





Transponder Depth by Cloverleaf 
(From Oceano [ 1985, p.12]) 


Maneuver 






Corrected 




Transponder 


Depth 



(m) 



1 


1366.9 


2 


1331.3 


3 


1 38’l . 5 


4 


1406.7 



C. BASELIHE CEOSSING METHODS 

Baseline length can be determined utilizing the baseline 
crossing method described by Haehnle (1 967) and Hart (1967) , 
which involves repeatedly transiting the baseline while 
simultaneously measuring ranges to two transponders- The 
minimum sum of the two horizontal ranges as determined from 
measured slant ranges and known transponder depths is taken 
to be the baseline length. 

The basic baseline crossing method requires a knowledge 
of: 

1- The ship’s course and speed, 

2- The accoustic slant range (or the two-way travel 
times from the ship to two adjoining ocean- 
bottom transponders, and 

3- The harmonic mean sound valocity- 
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The ship, travelling at a constant speed and heading, 
crosses the lines joining the transponders as shown in 
Figure 4.4. The acoustic slant ranges are recorded as func 
tions of time during the crossings. The sum of two hori- 
zontal ranges (the distance T1-T2 in Figure 4.4) is minimum 
when the ship is in the vertical plane containing the two 
transponders. Similarly, for a f our- transponde r array, the 
distances T1-T2, i'2-T3, T3-T4, T1-T4, T1-T3 and T2-T4 are 
determined. The minimum sums can be obtained by least- 
squares fitting the horizontal ranges against time and 
obtaining the vertex of the parabola as shown in Figure 4.5 
During the experiment four to six crossings were made for 
each baseline between transponders. The results are 
discussed below in Sections 1 and 2. 

Although the baseline crossing method is relatively 
simple with respect to data reduction, it is operationally 
time consuming. During the experiment it took one hour and 
thirty minutes for each baseline crossing. Thus the eight 
crossings (4 for north-south crossing, 4 for east-west 
crossing) took a total of 12 hours. 

1 - Msic Baseline Crossing Mo del 

Figure 4.4 shows a ship track crossing a baseline. 
The ship's transducer was set to interrogate the bottom 
acoustic transponders every 10 seconds for this experiment. 
Ranges were obtained from each transponder. It is easily 
seen (Figure 4.4) that, as the ship approaches the oaseline 
T1T2, the total horzontal range at the first ship point (0) 
is greater than at the second ship point (1) , and so on, 
until a minimum sum is reached at the instant of baseline 
crossing; from there on the total increases for each point 
on the track. Since the ship maintains a constant course 
and speed, the mathematical function which describes this 
relationship is a parabola (Figure 4.5) of the form (Hart, 
1967) i 
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Figure 4.5 Sample of Range-Time Parabola Showing 
Baseline Crossing Time. 

H = A T2 ♦ B T ♦ C (4.1) 

wnere 

A, B, C are constants to be determined, 

T rs ti me [ s ], and 

E is total horizontal range between transponders 
L m]. 

To determine the unknown variables A, 3 and C in 
Eguation (4.1) a least squares adjustment is performed. At 
least four sets of data R = R(T) are needed to solve 
Eguation (4.1). Since observed parameters contain random 
error. Equation (4.1) can be rewritten as: 



4 2 



A T2 ♦ B T + C 



H + V 



(4.2) 



where 



V is the residual of the observation. 



In matrix notation this corresponds to: 



(vr 






■ 


/ \ 










-t/'T..'/ 




A 

i 3 


= 


Rz 


; 




. 




c 




■ i 






^ ' • J 




j 







which may be written as: 



V + 
UK/ 



T 




3K/ 



R 



tlX.\ 



(4.3) 



(4.4) 



where 



V. 

L 




R. 

t 



are the residuals of the observations, 
are the numerical coefficients of the 
unknown variable, 

are the unknown variables (A, B, C) , and 
are constant terms of total horizontal range- 



The standard deviation of the mean of sound velocity 
is ±0.3 m/s (Wilson, 1960). So, the standard deviation of 
the distance is the following, where time is assumed error- 
less : 

d'R. = 0^.1 (4.5) 

Each observation is considered to be independent and 
uncorrelated. since a measurement of high precision has a 
small variance, and one of low precision has a large vari- 
ance, the higher the weight the higher is the precision and 
vice versa. Accordingly, the weight W of a single obsera- 
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tion is defined as a guantity that is inversely proportional 
to the variances of the observations (Mikhail, 1981, 

p. 6 6) . 



W = d2 / ds (4.6) 

where 6^ = variance of unit weight 

For the adjustment of Eguation (4.4), the least sguares 
criterion reguires minimization of the weighted function G 
(Mikhail, 1981, pp.6 9-73): 

3 = w, V,2 +W V2 ♦ ••• +F V2 = t^F.V.2 - (4.7) 

which in matrix form is written as 

i 

t 

G = V W V (4.8) 



where 

W = 



'\aJ,0 . ■ Q 
0 IaJ^ 



\ * 

Now Eguation 



ia/, 



rt 



nx/v 



(4.4) can be rearranged 



as : 



V = E - T X 



(4.9) 



Substituting for V into Eguation (4.8), get 

t t t t 

G = (R-TX) W(R-TX) = (R -X T ) H (R-TX) 



(4. 10) 



so. 



■f. t t t t t 

= R WR-X T WE-H hiTX + X T KTX 



(4. 11) 
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Since G is a scalar, the right hand side of Equation (4. 10) 
is also a scalar. Furthermore, the transpose of a scalar is 
equal to the scalar: 

= (x”^T^WR) = r"^WTX (4.12) 

and 



■t 

w = w 



(4. 13) 



Thus, Equation (4.10) can be written as: 



G = r'^WR - 2R^WTX + X^(t”^WT)X 



(4. 14) 



In Equation (4.14) all matrices are constants, except X, the 
matrix of unknowns. For G to ne a minimum, its partial 
derivative with respect to X must be equal to zero: 



^ t 

= -2R HT + 



9X 



2R"HT + 2X-^(T^HT) = 0 



(4. 15) 



and , 



(",.^ ,T ) X = T W R 
3xn tiM ^x3 3XJ 3X/t tun nx.1 



(4. 16) 



giving the solution for X; 



X = (T^ W T )-i W R 



(4.17) 



Now, the minimum distance occurs when dR/dT = 0 
which is given below: 



4 5 



T 



B / 2 A 



(4- 18) 



Thus, to get the aiaimum distance of the sum, one may 
substitute the time from Eguation (4. 18) into Eyuation 
(4. 1) : 



R (min.) = B2/4A - B2/2A + C 



(4. 19) 



where A, B, C are given by Equation (4.17) 



The variances of A, B and C (Mikhail, 1981, p258) 



are : 




(4. 20) 



The standara deviation of T in Equation (4.18) is: 



= [ (d^ B/2A2)2 ♦ (dg/2A)2]i/2 



(4.21) 



The standard deviation of the minimum distance K in Equation 
(4. 1 9 ) is;. 



2- Data Process ing 

During the experiment at sea, both north-soutn and 
east-west baseline crossings were run (Figures 4.6 and 4.7) 
using data collected by Oceano Instruments (1985a). 

For each baseline crossing, 11 consecutive interrogation 
points were fit to a parabola by least squares. The reason 
for using 11 points is that this number is sufficient to 
define the parabola near its minimum. The minimum value of the 
parabola can be considered to define the baseline length. 

The parabolas shown in the Figures 4.8, 4.9, 4.10, 4.11, 



d =( (T2d ) 2+ (id ) 2 +(52+ (2AT + B) 2(52 ]l/2 
R A o C T 



(4. 22) 
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4.12 and 4.13 (at the end of this section) are for single 
crossings between I1-T2, T2-T3, T3-T4, T1-T4, T1-T3 and 
T2-T4. Tne results are tabulated in detail in Table V, which 
gives the baseline lengths and their standard deviations for 
each crossing. To get the mean of these lengths, one can 
use the following equation: 



= (EH-&2^ ••• ^Bn) / n = f(E) (4.23) 

Since the standard deviation of the mean length E„in 

ffl 

Equation (4.23) is: 

6 ^ =[ (dl a;^/5E1)2+ (d2 a/^/aE2) 2+ ••• (4.24) 



4 7 



we nave 




Figare 4.7 East-iest Baseline Crossing. 



=[ 2+ (d2/n) 2+*.*+ (Cfa/n) 2 2 (4.25) 

Equations 4.23 and 4.25 are used with the data of 
Table V to compute the mean of baseline lengths and their 
standard deviations (Table VI)* 

Lyman BurJce of Oceano Instruments and L. Spielvogel 
of Seaco, Inc., have also computed the baseline lengths. 

They use a different method and technique using most of the 
data collected during the experiment. Their results and the 
differences between them are tabulated in Table VII- Due to 
our more limited data set and uncertainty in the depth for 
the baseline crossing method, these results differ from 
tnose obtained with the crossing method by about 2 m. 
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TABLE 7 

Baseline Lengths Determined bj Baseline Crossings (m) 



Crossing i 


T1 - T2 


1 


T2 - T3 


1 


T3 - T4 


1. 1 


1793.32±2. 43 


1 


1 993.09±0, 61 


1 


1972,39±2.50 


2- I 


1792. 62±2. 17 


1 


1 991. 92±2. 73 


1 


1969.50±2.02 


3. 1 


1793.67±2.27 


1 


1992. 75±2.38 


1 


1970.29±2, 19 


1 


1791.77±2.22 


1 


1 994. 33±2. 09 


1 


1969. 97±2.23 


Crossing i 


T1 - 14 


1 


T1 - T3 


i 


1 1 

1 1 

1 1 

1 Ki * 1 

1 1 

1 1 1 

1 

1 H3 
I ^ 

1 


1- 1 


1860.8 1±3. 25 


1 


2789. 62±2. 20 


1 


2588-66±2. 03 


2. j 


1858. 43±1. 81 


1 


2789.77±1. 62 


1 


2589.70±1.94 


3. 1 


1858. 63±2- 23 


1 


2789. 45±2. 1 1 


1 


2587.88±1. 14 


1 


18 58. 4 8±2. 08 


1 


2789. 73±1. 63 


1 


2588. 17+2.21 


5. 1 


1858. 65±2. 04 


i 


2790. 22±2. 26 


i 


2589.22±2.26 


6. 1 




1 




1 


2588-82±2.20 



3 . Mo difie d Baseline C rossi ng Method 

Mourad and Fubara (1972) developed this method. 

Their method determines the relative positions of ocean- 
bottom transponders and the depth of each transponder. 

There is no necessity to accurately measure each depth 
directly; this is the main difference from the previous 
work. The physical principle of this method is shown in 
Figure 4.14 S is the baseline length and Li and L2 are a 
pair cf coplanar ranges from the ship to transponders 1 and 
2 whose depths are Z1 and Z2, respectively. The mathmatical 
expression for this configuration is: 




riBLE ?I 

The Bean of The Baseline Lengths 





1 

1 


Mean Value 
( m ) 


1 Standard Deviation 

1 ( m ) 


T1 - 


T2 1 


1792. 84 


1 1. 14 


T2 - 


T3 1 


1992. 42 


I 0.93 


T3 - 


T4 j 


1970. 54 


1 1.12 


11 - 


T4 1 


1859. 00 


J 1.04 


T1 - 


T3 1 


2789. 76 


1 0.89 


T2 - 


T4 i 


2588. 74 


1 0.82 



S = ( L12 - Z12 )i/2 + { l 22 - Z22 )i/2 (4.26) 

To determine the three unknowns S, Z1 and Z2 by 
least squares method, at least 4 pairs of Li and L2 are 
needed. Now, Equation (4-26) can be rewritten as follows: 

F = S- ( L12 - Z12 ) 1/2 - ( l 22 - Z22 )i/2 (4.27) 

Assume 

A 1 = ( L12 - Z12 ) 1/2 (4. 28) 



A2 = ( L22 - Z22 ) 1/ 2 



(4, 29) 



Equation (4.27) is linearized and partially differentiated: 
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TABLE 7II 

Baseline Length Determination by Oceano 
Instruments and L. Spielvogel 



Lines 



Distance (m) 



Difference 



1 


Oceano 


1 Spielvogei | 


(E) 


T1 - T2 1 


179 2.68 


1 1791.7 


1 


0.98 


T2 - T3 1 


1997-46 


1 1995.9 


1 


1. 56 


T3 - T4 1 


1973.92 


1 1972-6 


i 


1.32 


T4 - T1 1 


1861-86 


j 1860.5 


1 


1.36 


T1 - T3 1 


2792-69 


1 2790.8 


1 


1. 89 


T2 - T4 1 


2592. 18 


1 2590.3 


1 


1. 88 



I 



J 



iL 

dU 



-A1- IL1 



9f_ 

QLZ 



-A2-IL2 



9F 

^S 

af 

aF 

9Z2 



1 



A1-»Z1 



A2-1 Z2 



(4.30) 



(4.31) 



(4.32) 



(4. 33) 



(4.34) 
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Figure 4.8 The Best Fit of Parabola for T1-T2 Crossing- 
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T2-T3 BASELINE CROSSING ! 





Figure 4.9 The Best Fit of Parabola for T2-T3 Crossing- 
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! T3-T4 BRSELINE CROSSING 




Figure 4. 10 The Best Fit of Parabola for T3-T4 Crossing- 
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TI-T4 BASELINE CROSSING 





Figure 4.11 The Best Fit of Parabola for T1-T4 Crossing, 
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T1-T3 BRSELINE CROSSING 
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Figure 4.12 The Best Fit of Parabola for T1-T3 Crossing. 
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2-T4 BASELINE CROSSING 
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Figure 4.13 The Best Fit of Parabola for T2-T4 Crossing. 
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/ j 

Figure 4.14 Baseline Crossing Configuration. 



Tne resultant observation equation for a leasr- squares solu- 
tion is (Mourad , 1972, p.21); 



AA + BV+W = 0 



{4.35) 



which can be written as: 



' / Ar'il, Aziz 






0 ■ 1 




Sj 


/ > 




" (9' 


« » • 


a5 




O . O -AU/j ■ 


Vi 




w'z 




0 


. 


A?) 

£,Z2 




. * * • 


'j 


+ 






• 


* • 4 

\ ) 


1 J 




1 • • . J 


IJ 








^0; 



Tnus, the xeas t- squa res solution (Mourad, 1972, 



(4. 36) 



p.21) gives: 
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(4.37) 



A = 







(BP;-‘b'^)-iw 



V =-P-ib'*^ (BP-iB^)-i (AA+W) (4.38) 

The variance of unit weight isc 

= {[ (BP-iB*') - 1 (AA+H) ]^W/df] 1/2 (4.3S) 

The variance-covariance matrix for the adjusted parameters Xq^ 
is: 



C(, £P^ +A^ (BP-ib'^^) - iA]-i (4.40) 

The adjusted parameters is: 

Xa = ^ 41) 

where 

A is a coefficient matrix for the unknowns correction, 

B is a coefficient matrix for the observation, 

V is a vector of residuals, representing the 
corrections to observed ranges, LI and L2 , 

A is the correction to assumed X^ , 

W is d constant matrix in Equation (4.27), when 
using assumed value (S<, ,L1 q,L2^,), 

P, is the weignt matrix of L^, 

P^ is the weighting function associated witn X^ , 
dg is tae variance of unit weight 
Xo. is true value of S, Z1 and Z2 , and 
Xg is an approximate value of X^.. 

The weight matrices are taken to be unity as the 
measurements were of equal precision. Table VIII shows the 
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resulting horizontal distances between transponders and the 
depths of each transponder. These results are obtained by 
measuring pairs of slant ranges at the instant the ship 
crosses the baseline. Since the ship crossed the baseline 
four times for the baseline crossing method, there are only 
four pairs of slant ranges which can be used; these repre- 
sent the smallest data sets needed for a least squares solu- 
tion. 

The baseline lengths computed by Oceano Instruments 
and Dr. Spielvogel are given in Table IX The difference 
between Oceano’s and Spielvogel's computed values are witnin 
the precision of the system. My results in Table 7III 
differ from TableIX, since Oceano and Spielvogel used much 
larger data sets (more than 250 data points) , whereas I used 
tne much smaller sets (11 points for each of four crossings) 
required by the two baseline crossing methods, namely eleven 
data points for each of four crossings. Since Oceano and 
Spielvogel used a resection method with four transponders, 
it is not necessary for the ship to travel the naseline 
repeatedly: The only requirement is that the ship oe in the 

vicinity of the center of the experimental area, where data 
can be collected continuously. On the contrary, the 
crossing methods require a large amount of ship time to 
collect a small amount of data. Thus, it is not surprising 
that there should be good agreement between Oceano and 
Spielvogel and rather poorer agreement between my results 
and theirs. 
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TABLE ?III 

Results of Hodifiel Baseline Crossing Method 



Trans- 

ponder 

pairs 


Horizontal | 

Distance Between | 
Transponders (m) 1 


Transponder 

Number' 


Mean 

Depths ( m) 


1-2 


1782.91±9.51 1 


1 


1361. 52±4.53 


2-3 


1993-02±6,30 | 


2 


1333. 94±2.54 


3-4 


1964.34±8.32 | 


3 


1375. 76±5.21 


1-4 

1- 3 

2- 4 


1836.65+9.21 j 

2815.47±15.60 1 

2574.13±7.55 1 


4 


140?’. 08±3.47 



1 










1 






TABLE IX 






Comparison 


of Besalts 


for Different 


Methods 










Basic 


Modified 




Lines 


Oceano 


Spieivogel 


Baseline 


Baseline 










Crossing 


Crossing 




T1-T2 


1792.68 


1791.7 


1792. 84±1. 14 


1782.91+9.51 




T2-T3 


1S97.46 


1995.9 


1992.42i0.93 


1993-02±6.30 




I3-T4 


1973.92 


1972.6 


1970-54±1. 12 


196 4. 34 + 8. 32 




T4-T 1 


1861.86 


1860.5 


1859.00±1. 04 


1836.65±9.21 




T1-T3 


2792.69 


2790.8 


2789. 76±0. 89 


281 5.47+ 15.6 




T2-T4 


2592.18 


2590.3 


2588.74+0. 82 


2574. 13+7.55 




1 










J 
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V. CONCLUSIONS AND RECOMMENDATIONS 



Two methods were used to determine the horizontal 
distances between the elements of an array of four acoustic 
transponders lying on the ocean bottom at an average depth 
of 1370 meters, namely basic and modified baseline crossing 
methods. The principal differance between the two tech- 
niques is that the former requires precise depth informa- 
tion, whereas the latter does not. A "cloverleaf" maneuver 
was described which may be used to determine the depths of 
the transponders. Because a larger number of acoustic range 
observations was used for the basic method, the standard 
deviations for the horizontal distances between transponders 
using that technique were about one-tenth those found using 
the modified method. Our results for the basic method are 
consistent with independent calculations made by Spielvogel, 
who used a much larger data set. The differences between 
our results using the modified method and Spielvogel 's 

are within our relatively large standard devi ati ons, wi th the 
exception of the distance measured between transponders 1 
and 3. This discrepancy is probably due to the very oblique 
crossings of the 1-3 baseline which were used. 

In addition to the smallness of the data sets used for 
our calculations other sources of error in our results may 
stem from the variability of the ship's speed and its pitch 
and roll, which we did not take into account. 

If the absolute positions of the bottom transponders are 
required, the position of the ship may be determined by an 
absolute positioning system such as GPS while measuring 
slant ranges from the ship to each transponder. When at 
least three or more absolute positions of the ship's 
acoustic transducer are known, the transponder coordinates 



63 



may then be adjusted by the method of least squares to 
obtain a best fit between measured and calculated slant 
ranges from the known transducer positions. This method is 
introduced in the Appendix. 
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AP PENDIX A 

ABSOLUTE POSITION OF OCEAN BOTTOM TBANSPONDEHS 



TLe absolute position of a ship can be obtained from the 
Global Positioning System (GPS). Here, ship position means 
the position of the GPS antenna. The antenna position is a 
function of time. Once the absolute position of the antenna 
is obtained, and simultaneously, the ship's speed, heading, 
pitch, roll and acoustic data, then the absolute position of 
antenna can be used to calculate the absolute position of 
the bottom transponders. The procedures are; 

1. lo determine the absolute position of the GPS 
antenna; 

2. To convert the absolute position of antenna to the 
hull-mounted transducer using the ship heading, pitch 
and roll data ; and 

3. Simultaneously, to convert the ship transducer's 
position to the nottom transponders by using acoustic 
data. 

A method xs introduced here to determine the absolute 
position of sea-floor transponders, if the absolute position 
of surface transducer is given. Least-squares adjustments 
must also be supplied to get the best fit ox the transponder 
position. 

The three-dimensional geometry used is illustrated in 

Figure A.1. There are 4 transponders (T1, T2, T3 and T4) . 

Their estimated coordinates are (X^,Y^,Z*), the absolute 

position of ship-mounted transducer is S(U;,V-,W. ), and tne 

J J J 

four slant ranges are 2*. 

Li 

The adjusted distances I* between transducer and tran- 
sponders can be computed given coordinates of any two 
points: 
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T4 




Figure A.l The Relationship of Transducer ( S ) 
To Transponders, ( ) . - 
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= [ <0-X1) 2+ ( V-I1) (H-Z1) 2 ]l/2 



(A. 1) 



5^ = [ (U-X2) 2+ (V-Y2) 2+ (W-Z2) 2 ]i/2 



(A. 2) 



= [ (U-X3) 2+ (V-X3) 2+ (W-Z3) 2 ]i/2 



(A. 3) 



% = [ (0-X4) 2+ (V-Y4) 2+ (H-Z4) 2 ]i/2 



(A. 4) 



which can be written simply as follows: 




(A.5)- 



wnere i are 1,2,3, 4, index corresponding to 
each transponder; 

j are 1,2 ■••• N. transducer positions 
at the surface; 

U^V, W are coordinate of transducer; 

X,Y, Z are coordinate of transponder; and 

are adjusted distance of slant range. 

It was necessary to determine the minimum number of 
positions (N) of the hull-mounted transducer that were 
needed to fix the position of four transponders in the 
three-dimensional coordinate system. There are three' 
unknown variables to be solved for each transducer position, 
and there are four transducer positions to solve for. To 
solve for these twelve unknown variables and apply least- 
squares adjustment, the number of equations must be greater 
than or equal to the number of unknown variables. Inis is 
given by: 
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4N > 12 



(A. 6) 



where 4 is four equations for each interrogation; 

N is number of interrogation; and 
12 is twelve unknown for 4 transponders coordinate. 

The solution of Equation (A. 6) is N>3- That is, the 
absolute positions in three-dimensional coordinates of all 
transponders can be determined with as few as three consecu- 
tive range measurements from four fixed transponders. 

Mikhail states the two-dimensional distance condition 
and its linearization (Mikhail, 1981, pp266-268) , which can 
be extented to three-dimensional coordinates. ■ Equation 
(A. 5) can be linearized: 



I.= io+[ (X.-U) /JO ]aX+[ (Y.-V) /J.o]ay+[ (Z.-W) /i.ojoZ; 

LO U t.t (.Li, C** 



The adjusted distance is then: 



J. = J. + V; 

C t U 



(A. 7) 



(A. 8) 



where J- are observed valuesof the distance; and 
are corresponding residuals. 

Equation (A. 7) then becomes; 



v-[ (X-U) /Jo ]AX.-[ (Y.-V) /Jo ]^Y.-[ (Z.-W) /J.o ]^Z = lp-S (a . 9) 

UC ct-t, Ltc ULI.'- 



Letting 

f. = Jo - (A. 10) 

in matrix form Equation (A. 9) becomes; 
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which may be concisely written as; 

V + B AX = f 
4/v/Xl 4N^X/2 /2X1 4^X1 



(A. 11) 



(A- 12) 



where 

V is residual of the observation; 

B is numerical coefficient of the unknown 
variable correction ; 

A X is unknown variable correction ; and 

f is numerical constant term . 

The weight matrix can be derived from the standard devi- 
ation of the slant range, which discussed in the previous 
chapter. The least squares solution of Equation (A. 12) is 
(Mikhail, 1981, p259): 

N = B^ KB (A. 13) 



t = B^ Wf 



(A- 14) 



AX = N-l t 



(A. 15) 



X = +AX 



(A. 16) 



V = f - BAX 



(A. 17) 
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(5 =[^2 (B N-» ) ]»/2 

A O 



(A. 18) 



^ t 

C2 = V H V / r 
0 



(A. 19) 



r = N - 3 



(A. 20) 



where 

N is coefficient matrix of the normal equations, 

W is weight matrix, 

t is vector of 'constants’ in the normal equations. 
Ax is vector of parameter corrections, 

Xois vector of approximate value, 

X is vector of adjusted value, 

d^is standard deviation of adjusted distance, 
o2 is estimate of the reference variance, 
r is number of statical degrees of freedom, 

N is number of given transducer position, and 

3 is number of observations necessary to specify 



Since during linearization we neglect all second and 
higher order terms, we must ensure that X is not signifi- 
cantly in error because of this approxima tion. 

Consequently , we must iterate the solution by using X as a 
new approximation and compute anotner correction- This 
procedure is repeated until the correction is insignif icant. 
The final result is: 



uniquely the model that underlies the adjustment 
problem. 




(A. 21) 



7 0 



wne re 



n is total iteration, 

is adjusted value at n iteration, 

Ax^ is correction at n iteration can be neglected compare 
with previous correction or within tolerance 
correction, and 

is true value of each variable. 

Since errors in transducer's position exist, it now 
must be established how this error affects the accuracy of 
transponder positions. This is discussed below: 

Let the standard deviations of transducer coordinates be 
r r Kecall Equation (A. 5) : 
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= [ (U 


-X. )2< 

c ' 


h (V 


-Y.) 


2 + 


(» -Z.)2]1/2 


(A. 


, 22) 
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unknown coordinate 


s.x,y 
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, are: 








= U 
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(v-yj 
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) 2 


]l/2 


(A. 


, 23) 


y. 
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= V 
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(W-Z. ) 


1 2- 


(U-X. 


) 2 


]1/2 


(A, 


, 24) 
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= w 


- [i.2- 
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(U-X. ) 

i. 


2- 


(v-y. 
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)2 


]l/2 


( A. 


25) 



Now assume 



Fl = [ J2.2- (V-y. ) 2- ( w-z. ) 2] 1/2 

L L i- 


(A, 26) 


F2 = [ $2- (W-Z. ) 2-( u-x. ) 2] 1/2 


(A. 27) 


c L ■ L 




F3 = [£2- ( 0 -X. ) 2-(V-y. ) 2] 1/2 

u C ^ 


(A. 28) 



So, the variances of Equations (A, 23), (A. 24) and (A-25) are 

(Mikhail, 1981, p.181): 

rfa =d2-(i. /FI) 2(<24- (v-y. /F 1)2^2+ (W-L. /FI) 2(j2 (A. 29) 

XUL g t ^ 
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62 =c[^2-(i. /F2J 2(^2+ /F2) 2(5^+ (U-Xj^/F2) 2d^2 



(A. 30) 



d2 =d^-(J?./F3)2(^2+ (U-X./F3)2(ii2+(V-Y. /F3)2d2 (A. 31) 

and the position accuracy is: 

= ± (d2 +d2 +d^2 ) 1/2 (A. 32) 
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